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ABSTRACT 

We study process of radiative neutrino pair emission \e) —>• \g)+ r y+vi> from coherently excited heavy ions 
(quantum mixture of two ionic states, the ground and an excited states) in circular motion. Determination 
of the neutrino mass is found to be possible with simultaneous detection of the photon and one of neutrinos 
in the pair down to the level of the smallest neutrino mass of order 5 meV in the three flavor scheme. 
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1 Introduction 


Conventional neutrino sources have been weak decay products of elementary particles (mainly pion and 
muon) and /3— nuclei. It was recently proposed [1] that circulating excited heavy ions may become another 
unique source producing large amounts of pairs of neutrino and anti-neutrino of the same flavor u a u a ,a = 
e, /j, t. Quantum mixture of excited and ground ionic states, realized by laser irradiation from counter- 
propagating directions against ions, inputs an ionic internal energy into the orbital motion, leading to a 
new kind of non-linear resonance. This resonance gives for rate calculations a stationary point in the crucial 
time integral of the phase factor, which otherwise monotonically varies, hence giving much smaller rates 
with much lower neutrino energies, as in the case of synchrotron radiation [2|. When the quantum number 
of excited level is appropriately chosen, neutrino pairs may be emitted with large rates during the resonance 
decay time. Produced neutrino energies may reach the GeV region suitable for oscillation experiments. 

In the previous work [3] physics related to the pair emission process |e) —>• | g) + vv was discussed. In the 
present work, we examine a different process: radiative emission of neutrino pair (RENP) |e) —>• \g) +7 + ^ 
from circulating excited heavy ions. Relative parities between the state \e) and the state | g) in this case are 
different (—), while they are the same (+) in the case of non-radiative neutrino pair emission. Thus, two 
processes cannot occur simultaneously. In order to avoid strong electric dipole transition (a purely QED 
process) against RENP, it is necessary to have a large angular momentum change like |AJ| > 2 between 
two ionic states, \e) and | g). A candidate ionic state for \e) is Ne-like heavy ion states in the configuration 
of (2p 5 3s)j=2, for instance Pb 72+ giving the first metastable state of J p = 2~ of energy 0(l)keV. In this 
work we shall not attempt a serious search for the best candidate ion, since this requires both of detailed 
experimental R and D works related to acceleration and dedicated atomic physics calculation. We shall 
concentrate on experimental principles and rate estimates assuming a hopefully reasonable set of parameter 
values. 

The rest of this paper is organized as follows. In the next section we recapitulate the main features of 
beam RENP from quantum coherent heavy ion beam, offering an intuitive understanding in terms of the 
non-linear resonance. Section 3 is the main part of the present work and provides calculation of RENP 
spectrum rates, both when neutrino variables are integrated out and when one of neutrino variables is left 
for its detection. Section 4 addresses the question of the sensitivity to the smallest neutrino mass in the 
three flavor scheme and discusses how small one can determine the smallest neutrino mass taking the case 
of normal hierarchical (NH) mass pattern.. 

Throughout this work we use the natural unit of fi = c = 1. 


2 Resonance condition of coherently excited ion in circular motion 

We first recapitulate the important feature of beam RENP. 

A promising way to prepare a coherent state of excited and ground states is via laser irradiation, often 
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more than two lasers for neutrino pair emission. It may lead to a quantum mixture of two states, 

|c) = cos 9 c \g ) + sin 9 c e lipc |e), (1) 

where the angle and the phase, 9 C and </? c , may vary in time, but we consider a shorter time scale than this 
variation time. Without any phase relaxation (a reasonable assumption after the circulating beam energy 
reaches its maximum), this state maintains its quantum nature, and various observable quantities are given in 
terms of the density matrix |c)(c| for the pure quantum state. Its off-diagonal element p eg = sm(29 c )e ~ llfc /2 
is called the coherence. An experimental method (adiabatic Raman excitation) realizing a high coherence 
is described and was successfully applied in [5] to gain a high level of quantum coherence of several % in a 
macroscopic system of hydrogen molecules in a gas phase. 

Unlike synchrotron radiation [2], the crucial time integral of phase factor econtains two contributions: 
the one coming from the ion circular motion and the other from the internal transition. These two terms 
may have phases of different signs in the time integral, and there is a possibility of phase cancellation. The 
cancellation is interpreted as a kind of resonance, and the resonance condition is given by the stationarity 
condition, dt& = 0. Emission of photon of energy uj and neutrino pair of energies E % , i = 1,2, in RENP 
gives the resonance condition, 

vt f 

(w + Ei + E 2 )( 1 - cos —) = — , (2) 

P 7 

when all three particles 7, u, v are assumed to be emitted at the forward direction. Since the angular 
distribution is well collimated around the forward direction, this is a reasonable assumption at least for an 
estimate. The constant 7 here is the boost factor determined by a circular velocity v: 7 = \fy/l — v 2 . The 
radius of the circular motion denoted here by p is equal to the Zeeman splitting energy eB/(M 7) under a 
magnetic field B for circular motion. The time t here is the look-back time measured from a fixed present 
time t = 0 , and for t -C p the resonance formation time t r is calculated from ([21) as 

t r ~pj-, -( 3 ) 

V 7(w + Ei+ E 2 ) 

The RENP photon spectrum is continuous and the total energy of three particles is bounded from above: 
u + Ei + E 2 < 2 e e 9 7 - This gives t r > p/ 7 . Numerically, taking the total energy of order its maximum 2 e e 5 7 , 
t r = p /7 ~ 3 x lCU 10 s for p = 4km and 7 = 5 x 10 4 . Dominant contribution to the phase integral comes 
from the resonance region, and this integral is approximately / 0 °° dt cos ((t — t r ) 2 /(2At 2 )). The width factor 
At r ~ y !p16 e g is of order 3 x lCU 12 s for e eg = IkeV, p = 4km, namely the resonance is very sharp, A t r <C t r , 
in the time domain, and is very broad in the energy domain. The resonance shape in the energy domain is 
not of a simple Lorentzian type. 

More precisely, the resonance behavior of the phase integral is worked out, using the technique of EL by 
the formula, 

[°° dt cos *(t), m = UJ + f + E2 ^D(t--D) 2 , D = / , (4) 

Jo 2 p 7 7 ui + Ei + E 2 

£ = — — ui — E\ — Eo + (1-^) 1 // 2 (o; cos ip cos 9 + pi cos ipi cos 9i) 

7 r i 
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for a neutrino pair emission of masses, mi,rrij. Emission angles (ip,0 for the photon, ipi,0i , i = 1,2 for the 
neutrino pair) are defined for the forward direction along the circulating beam to be at 8 = ijj = 0. For 
simplicity we took the small neutrino mass limit, which should be sufficient for our consideration of much 
larger neutrino energies. In this Airy-type of integral the development time t r and the resonance width Af r 
are 


resonance in time domain : 


t r = -D , 
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giving 


dt cos <3 >{t) 


width; Af r = m 


f°° , {t~t r ) 2 

/o (At r ) 2 


{to + Ei + E 2 )t 1 


> t r 


A t r = 


PI 


3 ^ to + E\ + E 2 


y./2 D -l/A 


( 7 ) 

( 8 ) 

(9) 


The allowed range of dimensionless function D is 0 ~ 1. The restriction D > 0 along with positivity of 
variables gives constraints on angular variables and energies of all three particles. 

It is necessary to irradiate lasers from the counter-propagating directions against the circulating beam 
at lease once in each revolution of circular motion, in order to efficiently re-excite lost available heavy ions. 
Moreover, irradiation from the counter-propagating directions effectively boosts laser frequencies in optical 
and infrared regions into the keV range of X-rays. 


3 RENP spectrum rates from circulating excited ions 

The RENP idea has originally been developed for small scale laboratory experiments aiming at systematic 
exploration of unknown neutrino properties [3]. The macro-coherent amplification (without the wavelength 
limitation) of otherwise tiny weak rates in atomic processes is the key concept for success of this project. 
This way one may effectively enhance weak rates by a factor oc n 2 where n is the number density of atoms 
in the available level. We recently observed the macro-coherent amplification in weak QED processes [5], 
but there are still many steps towards the final goal. 

In the present work we shall instead employ quantum coherence at a single ion level: the macro-coherence 
amplification is not necessary. But if it exists, it certainly helps much in giving much larger rates than given 
here. Acceleration of coherently excited ion in the present scheme is a new innovation and would require 
much R and D works, but it has a potential of producing coherent gamma ray beam if a macro-coherence 
is realized |TJ, which may have rewarding applications. 

We now work out RENP spectrum rates from circulating excited ions. The process for a single ion is 
l e ) - > Iff) + 7 + d = 1,2,3, where 1 7 and Uj are mass eigen-states (in the Majorana case Vi = Pj). 
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The important angular variable dependence of emitted neutrinos is in D , which is to a good approximation 
quadratic in angular variables as shown in eq. ([ 6 ]) . One can integrate over these angular variables by rescaling 
their magnitudes by I/ 7 , resulting in the following equation (11211 . The photon energy spectrum when two 
neutrino variables are integrated out is then calculated in a convenient integral form of two neutrino energies 
scaled by the maximum total energy co m , X{ = Ei/u m ,i = 1 , 2 , 

Um—r^- = RFij(-) , Uj m = 2 e e9 7 , (10) 

auj (j3 YYh 


R = 


2y/3(2n] 
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v 5 = f dV§{\ — r 2 ) 1/4 ~ 9.1x10 6 , Fij{y) = f dx x [ dx 2 H ij (y,x l ,x 2 ), 

J Jo Jo 

H i:j (y,xi,x 2 ) = y 5/2 ( 1 + ^-y)~ 2 xix 2 (xi + x 2 + y) 1/4 G ij (x 1 ,x 2 ,y) 9/4 Q(G ij (x 1 ,x 2 ,y )) 


„ , x , 1 ,m 2 rn 2 ; 

G ij (x 1 ,x 2 ,y) = 1-X! -x 2 -y- —^(-1--), 
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( 11 ) 

( 12 ) 

(13) 

(14) 


where | p) is an ionic state having a higher excitation energy than |e). In the example of Ne-like ion of 
\e) = ( 2 p 5 3 s)j = 2 one of these states may be | p) = (2p 5 3p)j = i. We assumed e pe = e eff /10 (close to Ne-like 
ion level spacing for this type of | p)) for simplicity. The dipole moment d pe is related to the El transition 
rate 7 pe which we assume of order 100 MHz. The step function 0(G) is defined by the property: = 1 for 
G > 0 and = 0 for G < 0. The neutrino energy range in the integral here is determined by G > 0 and 
Xi > mi/(2e eg 7 ). 

In this formula we incorporated finite neutrino mass effects only kinematically, hence there is no rate 
difference between the Major ana and Dirac neutrinos. The matrix element factor distinguishes Major ana 
and Dirac cases and shall be treated below. The boost factor ( 7 ) dependence is different from the case of 
neutrino-pair beam, and is larger by 7 2 than the previous case, as expected from a 7 — scaling law [T]. 

Without any experimental R and D works, it is difficult to estimate the number of excited heavy ions with 
a high coherence denoted by the parameter N\p eg (0)\ 2 . We took 10 8 for this value under this circumstance. 
This is the most important uncertainty in rate estimates of the present work. The overall RENP rate factor 
R of ea- dTIT) in the (e eg ,y) plane are plotted in Fig(Q[]). and Fig([2|) we illustrate photon spectral rates for a 
single massless neutrino pair emission, showing their sensitivity to the energy ratio e pe /e eg . We may judge 
from these figures that RENP from circulating coherent heavy ion can be measured with high enough rates, 
if N\p eg (0 )| 2 is not too small. 

A term 


3 mirrij 

8 EiEj 



3 mirrij 
8 oj% 1 xix 2 



TJ*.TT ■ 



(15) 


is added in the Majorana case to squared matrix elements common to Dirac and Majorana terms, c 2 ^ |, in 
the pair emission of mass sigenstates This term arises from anti-symmetrized wave functions of two 

neutrinos, only present for identical Majorana fermions [ 6 j. Here U e i,i = 1,2,3, is the 3x3 unitary matrix 
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Equi-rate, equi-energy curves 



Figure 1: RENP rate contours given by the overall factor R of eq. ffTTT) and equi-w m contours in (e eg , 7 ) 
plane. e pe /e e9 = 1/2. 7 pe = 100MHz assumed. 


Hz/w_m 


Photon spectrum 



Figure 2: Photon energy spectral rates for a massless neutrino: e eff /e pe = 10 in solid black, 8 in dashed red, 
and 5 in dash-dotted blue. Other assumed common parameters are N\p eg (0)\ 2 = 10 8 , 7 pe = 100MHz,p = 
4km, e eg = lkeV ,7 = 5 x 10 4 . 
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elements describing the neutrino mixing, and its elements have been determined from oscillation experiments 
[7] except one (in the Dirac case) or three (in the Majorana case) phase factors. The important difference 
from RENP of SPAN project |4J is that the conservation laws of energy and momentum do not hold, and 
contributions near the mass thresholds E) ~ m.j can exist for any given photon energy u. It is however found 
in the present work that the neutrino mass suppression given by rrii/Ei is severe at high neutrino energies 
Ej for the Majorana/Dirac distinction. 

Non-conservation of energy and momentum actually means that non-conserved amounts are compensated 
by ion de-excitation and its recoil, which are a small portion of their total amounts. It is however important 
to recover these by re-pumping ions with laser irradiation from counter-propagating direction at a straight 
section of the synchrotron machine. This way one justifies the method of rate calculation presented here. 

The main background when only a photon is detected is the three photon decay in which two photons 
escape detection. This background is rejected to some extent by subtracting escaped three photon events. 
The subtraction may be effective by different spectral shapes in the two processes. Still, there might be 
substantial background coming in after the subtraction. In this case one would have to rely on parity 
violating quantities such as emergent circular polarization that exist only in weak interaction process in 
order to unambiguously identity RENP. One needs a detailed study of the background rejection based on 
simulations taking into account a detector design. 


4 Neutrino mass determination and comparison with SPAN 


An efficient determination of the smallest neutrino mass m o requires a more refined measurement. It indeed 
becomes possible to determine the smallest neutrino mass, if one of pair neutrinos is simultaneously detected 
along with a photon. The double energy spectrum of 7 + 17 is given by 




d 2 V t 
' dudE h 


= RJ2 dx 2 H ij ( 


u E v 


U ri 


U r , 


,X 2 ) 


2 I r 3 m i m j «s ( : 
c ij \ - .2 „ : 


8ul,XlX 2 ' 


( 16 ) 


where 6m = 1 for the Majorana neutrino and 6m = 0 for the Dirac neutrino. 

Detection probability of a single neutrino event is estimated by the factor an^L where njv is the nucleon 
number density and L is the detector’s size along the neutrino beam. The cross section is of order 10 -39 ~ 
lCP 38 cm 2 for a 1 GeV neutrino, which gives cttinL ~ 10 -11 ~ ICE 10 (or smaller depending on lower detected 
neutrino energy) for a single neutrino detection using ~100 m detector size. One should multiply this factor 
to rate numbers shown in the presented figures, in order to derive actual detection rates of double events. 

To obtain realistic absolute rates for measurements, it is necessary to use a level spacing larger than of 
order keV, as is evident from rate numbers in Fig(JT]) and FigQ. Since we consider 0(l)keV energy range 
for the level spacing, e eg , the maximum neutrino energy 2 e e9 7 is at most of order 100 MeV for 7 < O(10 4 ). 
It is thus practical to consider double detection of combinations, 7 + u e and 7 + z/ e , since the charged current 
reactions Un —> ji require neutrino energies much larger than the muon mass. Thus, the relevant spectrum 
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rate is 


El^i 2 


d 2 Tj 
dujdE u 


(17) 


Neutrino detectors should be placed next to the photon detector at the synchrotron site to make easier 
coincidence experiments. 

The most promising case is to take the neutrino energy near the maximally allowed value = 2e eg r y, 
as shown for u m = lOOMeV in Fig([3|) and Fig(U]) in which the sensitivity to the smallest neutrino mass to of 
order 5 meV is indicated. The smallest neutrino mass range can be explored by devising a plot of deviation 
from the unity of the ratio relative to the three massless neutrino production rate, (c£T( 0 ) — dr(mo) )/dT( 0 ), 
which is shown for a few choices of the smallest neutrino mass m o in Fig([l|). The photon spectrum data 
may be obtained at different fixed neutrino energies, which make analysis of absolute mass determination 
easier. If calculated double detection rates are too small for a given detector design, one should think of a 
machine construction of larger boost factor 7 . An increase of 7 , for instance by 3, raises detection rates by 
a large factor, 700 in this example. There is no sensitivity to CP violating phases, unless two neutrinos are 
simultaneously detected. 


Double e rates: massive M 


Hz/w_m 



w/w_m 


Figure 3: Double spectral rates of photon with a detected v e of fixed energy 0.95c o m = 9.5MeV: with 
the smallest neutrino masses 5, 10,20, 50 meV’s all degenerate in the resolution of this figure. All cases 
for the Majorana NH of vanishing CPV paremeters,. N\p eg (0)\ 2 = 10 8 , 7 pe = 100MHz,/) = 4km and 
7 = 5 x 10 3 , e eg = IkeV, e pe = e eg /10 = O.lkeV are assumed. 


Let us compare the process of beam RENP discussed here with RENP of SPAN project [Jj. A marked 
difference in experimental methods is that the trigger laser irradiation is necessary for SPAN, while it is 
unnecessary for the beam RENP. Accordingly, the macro-coherence is required for SPAN, while it is only 
the coherence at the single ion level in the beam RENP. In terms of measurable quantities SPAN has a 
better sensitivity for the smallest neutrino mass determination because of many thresholds present, while 
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1-Double rate ratio: massive M 



Figure 4: Deviation of the rate ratio to the three zero mass case from unity for uj m = lOMeV: 5 meV in solid 
black, 10 meV in dashed red, 20 meV in dash-dotted blue, and 50 meV in dotted black, assuming the same 
set of parameters as in FigQ. The integrated deviation is 4 x 10~' 5 for the 5 meV case, while the double 
production rate in the threshold region of visible deviation is ~ 2.5 x 10 5 Hz. 


the threshold rise is hardly visible here. It would be instructive to write parameter dependence of rates in 
two cases. Disregarding common factors such as GpN\p eg (0 )\ 2 , the rate in the quantum beam scales with 
~ 10 -1 °7 6 7 pe (e eg /keV) 4 , while SPAN RENP scales with ~ l(n/10 21 cm -3 ) 2 7 p9 (eV/e eg ) 2 in the same unit. 
We took the radius factor of yjpteg = O(10 5 ) for this comparison. 

The important item for future R and D works is fabrication of (1) intense lasers with high quality to 
obtain large excited atom density n for SPAN, and (2) realizing heavy ion circulation with high coherence 
p eg for its RENP. In both cases a high coherence, either at the macroscopic level or at the single ion level, 
is of crucial importance for further developments. 

In summary, radiative neutrino pair emission from circulating excited heavy ions is useful for determina¬ 
tion of the smallest neutrino mass down to 5 meV level, if a high coherence can be achieved in the quantum 
heavy ion beam. RENP from circulating heavy ions is complementary to the neutrino pair beam, and two 
of them combined may give a comparative perspective to SPAN project. 
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process discussed in the present work. This research was partially supported by Grant-in-Aid for Scientific 
Research on Innovative Areas ’’Extreme quantum world opened up by atoms” (21104002) from the Ministry 
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